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III! ('ONTRIUUTION Ol- Till-; Dll-rUSli LKiMT COMPONENT 


1C) Tin: TOPocnuPiiK' i-i fi-ct on rf:moti:ly si-:nsf:i) data 

Chris Justice 

Brent llolben 

ABSTRACT 

The topographic ctTect is measured by the difference between the global radiance from inclined 
surfaces as a function of their orientation relative to the sensor position and light source. The short- 
wave radiant energy incident on a surface is composed of direct sunlight, scattered skylight and ligirt 
refiecteil from surrounding terrain. Tlte latter two components arc commonly known as the diffuse 
component. I'he objective of this study was to examine the contribution of the diffuse light com- 
ponent to the topographic effect and to assess the significance of this diffuse component with re- 
spect to two direct radiance models, to spectral band ratioing and to simulated Landsat data. 


Diffuse find global spectral radiances were measured for a .series of slopes and aspects of a uni- 
form sand surface in the red and photographic infrared parts of the spectrum, using a nadir pointing 
two-channel handheld radiometer, The diffuse light was found to produce a topographic effect which 
varied from tlie topographic effect for direct light. The topographic effect caused by diffuse light was 
found to increase slightly with solar elevation and wavelength for the channels examined. The correla- 
tions between data derived from two simple direct radiance simulation models and the field data were 
not significantly affected when the diffuse component was removed from the radiances. Diffuse rad- 
iances contributed largely to the variation in ratioed data. Subtraction of the diffuse radiance prior 
to ratioing resulted in a 50 percent average decrease in the standard deviation of the ratioed data. 


In an extreme case of radiances from a 60 percent retlective surface, assuming no atmospheric 
path radiance, the diffuse light topographic effect contributed a maximum range of 3 pixel values 
in simulated Landsat data from all aspects with slopes up to 30 degrees. Such a variation is suffi- 
ciently small compared with other variations which arc likely to occur in the data, to indicate tha' 
the diffuse component docs not warrant modeling prior to cover classification analysis. 


rm- (ONTRiBUTioN ()i riU' Dii'i-usi-: ur.iiT componi-nt 
TO Till TOPOCRAFIIK' I Fl'IXT ON RliMOTI-LY SI-INSI I) DATA 


1. INTRODUCTION 

The shortwave radiant energy inciilent on a horizontal portion of the earth’s surlaee is eomposed 
of direet sunlight, scatteretl skyliglit and liglit reileeted from surrounding terrain, The hitter two 
eoni|>onents are in common usage termed the diffuse component (Lui and Jordan, I960; Stanhill 
1966: Kondratyev, 1977), It is impoitant to note that certain studies have used the term diffuse 
light to refer strictly to unpolarized light (Shurcliff, 1962), I hroughout this study the common 
usage was adopted. 

In modeling the sensor response from inclined surfaces, llolben and Justice (1979), Justice and 
Ilolben ( 1 979), and llolben and Justice ( 1 980) showed that there was a need to examine the contri- 
bution of the diffuse light radiance to the topographic effect. The topographic effect is defined here 
as the variation in global railiance from inclined surfaces as a function of their orientation relative to 
the sensor position and light sources. The term “topographic effect” is used because in the context 
of satellite remotely sensed data e.g,, Landsat data, the light source and sensor geomet' v are e.ssenti- 
ally constant for a specific data set and the principal variation is in the surface geometry, i,e, topog- 
raphy, The degree of topographic effect can be measured by comparing the radiance from a hori- 
zontal surface to the radiance from a sloping surface of the same cover type. The objective of this 
paper was to quantify the contribution of the diffuse light component to the topographic effect 
through field measurement and to assess the significance of the diffuse light contribution both for 
two selected theoretical models and simulated Landsat data. 


2, BACKGROUND 

Tire irradiance received at a surface is the sum of the direct and diffuse components and is often 
termed the global irradiance. The diffuse skylight component impinging a surface is caused by the 
atmospheric Rayleigh and Mie scattering of the solar beam and varies in proportion to the direct 


conipoiuMU as a function of wavelength, optieal path length, eoinposition of the atmosphere (Bull” 
rich I ‘)(i4. anil l■l•asel• I ‘>75) aiul the orientation of the surfaee. Measurements of the intensity of 
iliffusc radiation ineiilent on a hori/onlal plane at the l arth’s surface and the proportional relation* 
ship to the global nu.\ have been presented both for elear skies (h]bal I ‘)70a and b) and cloudy con- 
ditions (.Stanhill I %.'>). Measurement of the diffuse skylight is usually taken by obseiirring the solar 
disc, although the effeetiveness of this method is questioned by Ileyvvood ( l‘)6b). Measurements of 
the spectral distribution of both direct and diffuse components at varying times of year are presented 
by Boer ( 1U77), Hourly, daily, and monthly totals of diffuse radiation calculated from surfaee 
measurements were described by Lui and Jordan ( l‘)bO), Norris ( l%b) and Ooldberg et al ( l‘)7d). 
Diffuse and direct irradiance measurement have also been made for sloping surfaces by using hemi- 
.spherical pyronometers tilted at various angles and aspects (Kmulratyev and Manolova, l‘)6(); lley- 
wood, I ^65; Temps and ('oulson 1 d77; Klueher 1 47‘)). By definiiion, the diffuse skylight incident 
on a surface is multidirectional. The path length of atmospheric scatterers such as aerosols and air 
molecules changes with viewing direction and therefore tlie spectral distribution of skylight intensity 
would not be expected to be equal around the celestial hemisphere, This anisotropic distribution 
is important when consiilering the proportion of the dlffii.se component relative to the global ir- 
radiance impinging inclined surfaces, 

1'he general distribution of tlie inten.sity of sky radiation has been presented by Bullrich (1064), 
Bullrich et al. ( 1 068), Coulson (1071) and Temps and Coul.son ( 1 07 7) and shows an ani.sotropic dis- 
tribution around the celestial hemisphere with a maximum around the solar aureole, caused by strong 
forward scattering and a minimum, normal to the solar beam in the principal plane, There is also a 
marked increa,se in brightness towards the horizon, due to limb brightening of the l‘art'i, Temps 
and Coulson ( 1 077) observed that for a 34° solar elevation that the intensity of skylight was 40 per- 
cent greater at the horizon tliat at the zenith. 


I’edorova ( 1 065) (referred to in Kondratyev 1 077) showed that 80 perc 'Ut of the total scat- 
tered flux incident on a horizontal surface will come from the circumsolaf iialf of a clear sky. The 
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ink'iisity ut'ilitTiisc lijihl reaches a peak hetween .4 am! ,5^m due to molecular Kayleiph scattering 
((lates 1%5. Dave etai. 1 475) and the ratio of the irradiances ol‘ direct sunlight and global llux decrease 
with decreasing wavelength (Dave etal. 1475). The dilTuse component ofirradiance incident on a liori- 
/ontal surface is greater than the ilirect component for wavelengths shorter than JlSixnx. I’he intensity 
and distribution of skyliglit aroumi tlie celestial hemisphere is sliown to vary lH>th with solar elevation 
and wavelength (Volt/, and Uullrich. |46| , Dave 1478, Hullrich et al. 1478). 

Biillrich et al. ( 14o8) sliow that tlie skylight minimum becomes angularly more distant from the 
sun with increasing solar elevation anti wavelength. They reporteil that the skylight radiance iliffer« 
ential between the solar /enith and the horizon increases with wavelength. 

The relative proportion of iliffuse to global irradiances increases ilirectly with solar zenith angle 
(b’raz.er 1 475), this increase being a maximum at large incidence angles (Coulson 1 47 1 ), Temps and 
(’oulson (1477) showed that under midlatitude, contenental summer, clear sky conditions, the sky- 
light contributed 1 6 percent to the total flux on a horizontal surface, for a solar elevation of 37 tle- 
grees, Hey wood ( 1 466) states that the diffuse component on a clear sky midsummer day will rise 
from 16 percent of the total radiation at noon to 25 percent at 5 hours before or after noon. Dave 
et al, (1475) stated that the diffuse component on a horizontal surface is about one seventh of the 
ilirect component for the sun overhead but 1 .2 times as much as the direct cominment for a solar 
zenith angle of 80°, under average midlatitude summer conditions. 

Most studies of diffuse radiation are undertaken under clear sky conditions with no apparent 
visible haze. A change in aerosol content, water vapor and/or cloud cover will however cause a 
change in spectral intensity and spectral distribution of the skylight and a change in the relative pro- 
portions of direct and diffuse light. Detailed discussions of the intensity distribution of the .skylight 
under different atmospheric conditions are given by (’oidson (1471 ), Dave ( i 478) and Dave ( 1 474), 

The component of the diffuse light scattered from adjacent terrain contributes to the global 
llux incident on a surface but has received little attention in the literature due to the complexity of 
isolating this component for measurement. Light renected from adjacent surfaces either onto the 
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ineasureU surlace. or into the sensor Cielil of view will be dependent on a lunnber of factors, such as 
directional rellectahce characteristics, albedo and orientation of the adjacent surfaces. Ileywood 
(l‘)00) stateil that terrain rellection and reradiation is neglijiibie on surfaces of less than 60° slope, 
but can attain about H percent of the global radiation on a vertical surface in the summer. 

Boer ( 1 077) gave a general estimate of the contribution to the total sunlight by renectance 
from ailjacent surfaces of between 4 aiul 7 percent for surfaces normal to the solar beam, assuming 
a 20 percent surface albedo. Kondratyev ( 1 077) staled that for steep slopes with high albedo values, 
the reflected radiation from surrounding terrain may constitute a considerable proportion of the 
global llux and that tliis proportion will be iiighest at low solar elevations for slopes lacing away from 
the sun, In an example for a surface with a 20 percent albedo and in an extreme case (i.e,, low solar 
elevation for slopes away from the sun), he reports that the terrain renectance can constitute 6‘) 
percent of the global Hux for a ‘H)° slope and 9,8 percent for a 30° slope. 

Modeling the radiation incident on surfaces has been undertaken for a number of applications. 
Recent developments in modeling the radiation from surfaces have been made for remote sensing 
applications (e.g., Oliver and Smith, 1 974, Marks and Dozier, 1 979, Kimes and Kirchner 1 980), 
whereas previous insolation modeling had predominantly meteorological and solar energy applica- 
tions. Several of the early models assumed an isotropic distribution for diffuse light under clear sky 
conditions. Lui and Jordan (1 963) developed an insolation model to predict the global radiation 
on inclined surfaces assuming an isotropic diffuse sky distribution, llolben (1 975) developed a 
direct and isotropic sky diffuse insolation model which incorporated first order effects due to slope 
orientation and shading from adjacent topography, However, several studies have shown the adverse 
effects of applying the isotropic diffuse sky assumption (e.g., Temps and Coulson 1977; Kluchcr 
1 979; Dave 1 979). Kondratyev (1 977) shows convincingly that the isotropic assumption is 
accurate for slopes facing perpendicular to the principal plane and for slopes of le.ss than 45°, 

The general anisotropic “all-sky” model presented by Klucher (1 979) provides a good model 
for thesolar radiation on tilted surfaces. The model was developed from horizontal surface measure- 
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HK'iUs aiul was sliovvn to liavc a systematic error oDv-ss than 2.5mw/em“. I'liis moild showed im- 
provements over the models presenteil hy Liii and Jord-”' ( 1 ) and Temps and Coulson ( l‘)77). 

Inclusion of the scattered terrain radiance in insolation nuulels <)ii inciineil surfaces has been under- 
taken hy Kondratyev ( l‘)77) Dozier ( l‘)78) and Kiines and Kiivhner ( l‘)H0) hut these models are 
suhstantially more complex, Dozier ( lh7S) developed a comi’rehensive .solar railiation model for 
snow surfaces In mountainous terrain and in a .sample of the model output of a typical mountainous 
situation, showed the following results: for the ultraviolet and visible spectrum 75 percent of the 
incoming radiation was ilirect. 15 percent iliffuse anil 10 percent from terrain renectance: for the 
infrared h| percent was direct, 2 percent diffuse, and 7 i»ercent terrain renectance. Kimes and 
Kirchner ( l‘)80) demonstrated that exclusion of adjacent terrain renectance when modeling 
Lambertian target rellectance under extreme conditions, i.e. adjacent slopes of 30°: sun angle cor- 
rected albedo of 5 percent: infrared wavelength, could lead to a d percent error, 7*he error tc<‘in was 
considerably less for slopes undV- 30°. 


For remote sensing studies and in particular examination of the topographic effect sv are 
concerned with radiance emanating from the surface and not solely the insolation or irradiance im- 
pinging the surface and thus we need to consider the reflection of both the direct and diffuse com- 
ponents. The radiance measured from a surface is a function of the reflectance properties of the 
surface, the atmosphere between the ground and the.sen.sor. the .sensor geometry, and the incident 
radiation. Fxamination of the diffuse radiance component in the context of remote sensing will 
help us to model the topographic effect more accurately and to assess its significance in models ap- 
plied to satellite remotely sensed data, F’or non-Lambertian .surfaces, the radiances will be determined 
by the directional scattering properties of the surface, the intensity and direction of radiation 
sources, and the exitance angle and therefore radiance modeling will undoubtedly be more complex 
than for Lambertian surfaces. 
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(iROUNDDA I ACOI U cnoN 

TIk* iipi’roacli uiloptod lor quantilVinp the ililTiisv light component to tlie topographic elTect 
was to measure ililTu,-. ami global radiances from surfaces tilted at a range i)f slopes ami asix'cts. 

’a 'iances were measured consecutively fumi a imn-Lambertian coarse samI surface and a BaSO^ 
reference surface. The BaS ()4 surface approximates a 1 amhertian surface and has a rellectance of 
approxiniately d8'< . for the wavelengths and view angles examined in this study. At ste.’pcr view 
angles a BaS().| surface approximates a Lambertian surface far less (Schutt |d7(>, llsia and Richmond 
|d7o). The I ambertian reference plate provides a measure of the irradiance on the surface and can 
be used to calculate the surface reflectances. 1'he reflectances calculated for different light source 
ami sensor geometries are termed reflectance facttirs (Judd, l‘)o7). I'he term bidirectional renectance 
factor is more commonly used (Kriebel, I d77. Nicodemus et al., I d77) and for a specified wave- 
length gives a measure of the amount of light reHected into the sensor, relative to the total amount 
of light impinging the surface. More detailed discussions of the bidirectional reflectance factors from 
natural surfaces are given by ('oulson et ah. ( l%.‘') , Coulson ( ld()()), Kriebel ( lh7(), 1 077. 1 078). 
and Robinson and Biehl ( 1 070). 


Radiance measurements were obtained for both .surfaces, inclined at all combinations of slope 
angles, ranging from 0-60 degrees in lO” increments and aspects for the 1 6 compass points, in 22 5^’ 
increments. The surface aspect was measured in degrees clockwise from the sun’s a/.iimith. This 
angle is termed the "azpect" of the surface (llolben and Justice, 1070). Red am! I’hotograpliic in- 
frared radiance data pairs were collected in data subsets called “a/pect strings,” that is slopes be- 
tween 0 - 60° in 10° increments for each azpect. The measurements, in ir radiances, were taken using 
a two-channel nadir pointing handheld radiometer, similar to that described by Pearson et al. (1070). 
Illtercd for the red (0.63 - 0.60 jum) and photographic infrared (0,76 - 0.00 nm) bands. These chan- 
nels are cciuivalent to the proposed 'fhemiitic Mapper bands 3 ami 4 of Landsat 0 ('fueker et al., 
1080), Although n radiances were measured using the handheld radiometer, the term radiance is 
used to describe these ineasuremenirthroughout the text. 
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riic surfaces surrouiuling tlie target were painted tlat black to minimize terrain scattering. (»Io* 


bal and diffuse radiances were measured for both surfaces at each slope/azpeet combination. Tlie 
diffuse observations were obtained by shielding the solar disc with a small opaque panel. Three data 
sets were taken under clear sky conditions for solar elevations of 23“, 2^°, and 3^5'’. which correspond 
to the Laiulsat sensing time for midlatitude solar elevations occuring during late fall, winter ami early 
spring, or for high latitude solar elevations occuring iluring late spring, summer and early fall respec- 
tively. I^ach data set was collected in less than 45 minutes and each azpect string in less than 3 mln« 
utes. The measurement apparatus was reoriented to the sun's azimuth after collection of ratiiance 
measurements for each azpect string subset, to reiluce errors due to the apparent movement of the 
sun. 


4. ANALYSIS 

The analysis of the data is presented in four subsections. The first section (4.1 ) describes the 
general characteristics of the diffuse light topographic effect. Section 4.2 shows how the diffuse 
light topographic effect changes with solar elevation and wavelength. Section 4 ■ ;sess'*s the sig> 
nificance of the diffuse component with respect to two direct radiance moilels and to spectral band 
ratioing. The final section (4.4) demonstrates by a simulation study, the contribution of t/n diffuse 
light topographic effect to Landsat MSS sensor response. 


4. 1 General Characteristics of the Diffuse Light Topographic Ld’fect. 

A topographic effect is observable when a change in radiance occurs due to a change in surface 
orientation. This topographic effect can be (luantified by .subtracting the radiance measured from a 
horizontal surface from the radiance measured from an inclined surface. Variations in diffuse radi« 
ance from tiic BaS 04 surface were observed between slope angles and orientations (L'igure 1) and 
therefore by definition a topograpliic effect was apparent due to the diffuse component, The great- 
est radiance range occurred in the principal plane of the sun and the smalles'. range for slopes perpen- 
dicular to the principal plane. An increase in radiance with slope occurred for all azpects. This pat- 
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Icni strongly resembles tlie anisotror-ic ilistribution of ililTuse liglit around the celestial hemisphere, 
descrilK'd in Section 2. Tlie variation in the diffuse radiance values is directly related to variations 
in the sky brightness, the maximum sky ratlianee in the .solar half of the celestial hemisphere ex- 
plains tlie greater radiance values associated with a given slope angle, for those azpects close to solar 
azimuth. 7'he strong positive relationship between slope ami radiance can be explained by the limb 
brightening of the Harth and a possible contribution from terrain reflectance at the higher .slope 
angles. The close correspondence between the anisotropic sky distribution and radiance would be 
expected iTom a Lambertian surface, as the radiances are directly related to the intensity of the im- 
pinging radiation and are independent of the view angle. 


All natural surfaces are non-Lambertian and to some degree have preferred orientations of 
scattering, i.e. angular anisotropy of renection (Kriebel IO 76 ), The radiances from non-Lambertian 
surfaces are therefore dependent on view angle. The degree of departure from Lainbertianess deter- 
mines the magnitude with which radiances will be affected by view angle. The diffuse radiance values 
for the coarse sand surface in tlie plioto, raphic infrared (0.76 - 0.90 ^m) (Figure 2), reveal a .similar 
pattern to that displayed by the Lambertian radiances (Figure 1), namely a marked increase in radi- 
ance witli slope, Kxamination of the radiances for the same slope angle with different azpects (i.e., 
tlie concentric circles in Figure 2) reveals little variation, indicating that the scattering from this non- 
Lambertian surface is not particularly orientation dependent. This is most likely due to the random 
distribution of the sand grains and tlieir individual reflecting surfaces. The similarity between the 
coarse sand and BaS 04 surfaces indicates that the sand surface has no extreme preferred direction(s) 
of scattering. This is confirmed by examination of the direct and diffuse light reflectances for each 
slope and azpect combination (Figures 3a and b), The reflectance as shown in Figures 3a and b is 
die ratio of the radiance from tlie non-Lambertian surface (coarse sand) to the radiance from the 
Lambertian reference surface (BaS 04 ) for a given slope and azpect. The relatively constant reflec- 
tance values for the same slope angle at different azpects (Figures 3a and b) show the lack of orien- 
tation dependence of scattering from the sand surface. The difference between the reflectances in 
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I'ijum’ .Ml iiulicali.’s (luit Iho isotropic assumption l\'r tlto ililTusc liglit will not accurately model the 
radiance. .Mthounh in this case the natural suil’ace examined did not show any i tarked directional 
scattering, it can be hypothesized that the variations in sky brit’htness combined with the reflectance 
properties of the surface will make diffuse radiances from natural surfaces cmnplex to nurdel. 


d. ’ N'ariations in the Diffuse I ipht I'opopraphic I ffect with Wavelength and Sidar 1 levation. 
rite diffuse lii’ht topoi’.raphic effect for the Ha.SO,, surface in the red part of the spectrum 
- 0.(P) jam) illustrates the same I’.eneral pattern in the radiances (Fipure 4) as for the infrared 
(|‘i!',ure 1 ) portiim of the spectrum. Both the diffuse radiance values and the ranp.e for each a/pect 
were lower for the reil than the photographic infrarcil, the wider pliotopraphic infrared bandwidui 
examined in this study, resulted in hipher overall radiance values. The dccrca.se in the radiance values 
for 0-20" slopes facinp away from the sirlar a/imulh. is cau.sed by the dominance of the skylipht 
minimum radiance in the cxpo.sed portion of the celestial hemisphere, rel itive to the coirtribution 
from the hori/on and surroundinr. surface reflectance. 


Variations in the Ba.SO, diffuse lip.ht radiances were pre.sent for all data .sets collected, ’ amples 
for .’.I" and .solar elevations in the red and infrared wavelengths are presented in l■i!Ulres .‘?a. b, 
c, and il. In p.eneral, hip.her solar elevation data sets had hip, her diffu.se radiances and sliphtly preater 
ranpes in the radiances associated with an a/pect strinp. particularly in the photopraphic infrared 
( fable 1). fable 1 shows a slipht increase in the diffuse lipht topopraphic effect, fhe increase in 
the radiance with solar elevation is ilue to the smaller anple of incidence made with the surface 
and the resultinp increase in the plobal flux throuph a decrease in atmospheric path lenp.th. It should 
be noted that althouph the diffu.se radiances increase with solar elevation, the proportion of the 
diffuse contponent relative to the direct decreases, 


(’alculated differences in the toi'opniphic effect for three solar elevation data sets, as measureil 
by the ranpe in radiances associated with each a/pect, showed that the preatest variation in the raili- 
ance ranpes between data sets, occurred for those slopes facinp into solar a/imuth ( fable 1 ). fhe 
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loiHH'.iai'liic olTool was ossculially IIk> samo for ilu>sc a/|'oois I.K'liSf’ away iVom solar a/iimHh ami 
vaiicvl loss ri>r (ho loil ohaniiol. 


•}..? I ho Kolalivo IVoportion o{' (ho nuTuso I iphl ('oiuponoiK ti> (ho Global Kavlianoo 
Koioivoil by (ho Sousvm. 

In (ho provious (wo soo(ions, (ho analysis iloall wi(h Iho inapaiiuulo ami variabilily o( (ho ilib 
Cnso I4',h( lopof',ia|>hio olTool I'ui (o nmlorsiand (ho oondibniion oC (his oomponoiK (o (ho lopoiiiai'hio 
olTool, i( is nooossary (o oonsiiloi (ho maj’nilmlo oC (ho iliCCnso lailianoo rolalivo lo (ho global railianoo. 
Iho amoniU o( MilTuso lighi loooivoil by (ho sonsor rolalivo (o (ho global radianoo (Voin ll»o UaS()_, 
snrCaoo was oalonlalod as a poroonlago ami is plodoil by slopo ami a/pool foi (wo solar olovalions, .d>" 
anil 'X' (I'ignros oa ami ob rospoo(ivoly), I illlo ohango wasobsorvoil in (ho proporlion oCililTiiso 
liglK roaohing (ho sairl'aoo Cor (hoso a/poois I'aoing in(o solar a/imn(h, allliongh a niarkoil inoroaso in 
Iho proporlion of (ho ililTnso ooni|H>non( ooonrroil wi(l( sloi'o, Cor (lioso a/pools porpomlioiilar lo ami 
awa^ iVoin solar a/imuih {I'igiiro oa ami nb). A inarkoil inoroaso ooonrs as inoiilonoo anglos of 
ai'o ai’proaolioil (i.o.. a( grazing anglos). I'or a/|'oo(s I'aoing in(o sun (boro is a gonoral ilooroaso in 
I’oroonlago ilii'Cnso lighi wilh incroasing slopo. This is ooimloraoloil by (orrain rol'loolanoo a( high 
.slopos, i.o. groalor (han 40". I'ho samo gonoral I'aKorn oan bo .soon Cor (ho (wo solar olovalions U'ig- 
ni'os na anil ob), Inil il shoulil bo noloil (ha( inoro slopos aro in shallow Cor (ho ilala sol, Tho ilala 
Cor lowor sim olovalions .showoil a smallor porooni iliCCn.so oompononi Cor slopos I’aoing inlo sun ami a 
largor porooni Cor (hoso I'aoing away Crom solar a/inuilh. I'ho proporlion oC iliCCu.so lighi Cor (ho hori” 
zonlal snrCaoo is similar Cor bolh ilala .sols (i.o., o. 1.^ porooni). I'ho varialion in railianoos wilh azi'ool 
Cor (ho hori/onlal snrCaoo is oau.sod by slighi ohangos in solar olovalion ihiring (ho moasnromoni pro- 
ooiluro. I'ho (opographio oCCool oansoil by (ho illCCnso oompononi was moro pronoimooil a( low .sun 
anglos, (han a( high sun anglos (b'iguros (>a ami ob). .Allhough (ho absolulo iliCCu.so railianoos aro 
smallor al low sun anglos, (ho proporlion oC (ho lolal inooming iliCCu.so railialion is gonorally highor. 
iho ooniribulion oClho iliCCuso oompononi (o Iho (opographio oCCool varios oonsiilorably wilh aziH'ol 
(I'iguros Oa ami ob). .AKhough Iho iliCCu.so railianoos vary rolalivoly lilllo oomparoil (o (ho ilirool 


liuliiliu’v's, llu' ivlativo proportions of each componont vary oonsiiloraNy. iiulicating that an isotropic 
assumption lor ililYiisc liphl svonlil fit llic ticUl ilala with ililTorin{i ik'ijcos of success, ilcpciuling on 
slope aiul a/pcct conripiiralion aiul solar elevation. 


4,4 .Vssessini', the Sipnificanee oT the DilTiise t'oinponeni lor I'wo Direct Kailiance Moilels 
aiul Spectral Ikiiul Ratioinp. 

The sipniHeance of the ililTnse compos cm only he a.scessetl with reference to a piven applica- 
tion of the railiance ilala. In this simly, we aic v.siuvmevl with the sipnificance of the iliffnse compo- 
nent in moilelinp railiance data to eliminate the topopraphic effeel. One way to evaluate this sipnifi- 
canee is to examine the effect of the iliffnse component on two previously nseil ilirect railianee moil- 
els (.Instice aiul llolhen, hy correlalinp Iheoretieal railiances from the railiance moilels to the 

measiireil plohal railiances and the calculated direct radiances, .V second assessment of the sipnificance 
of diffu.se radiation is undertaken hy examininp the effect of the diffuse component on spectral hand 
ralioinp. 


I'he sipnificance of the diffuse component on the correlation helween the field measured data 
and simulated data derived from sunliphi models proposed hy .lusiiee and llolhen ( l‘)7‘)) was a.s- 
sessed. Data derived from the simple 1 amherlian model (C'os i) was correlated with the p.lohal radi- 
ance for the HaSO,, surfaee. fora sun elevation of .to", piu. direct radiance was calculated hy suh- 
Iraclinp the diffuse radiance from the plohal radiance, .A small improvement in the coefficients of 
determination (r’) of up to o pereent was obtained for the radiances with the diffuse component 
suhtracled. I'he improvement was partieularly marked for those a/peels perpendicular to solar 
a/inuith. ('orrelation coefficients for the data derived usinp, a non-l amherlian model, ('osl'i • Cos*' ' e 
(.lusiiee and llolhen D)7‘). Smith el al. l‘>80). were caleulated with the plohal and direct radiance. 

An averape improvement of 7(i percent in the eoeffieient of determination for the radiances with 
the diffuse component suhlracted over the coefficients for the plohal radiance was ohserved for those 
a/pecis iterpendicular to the .solar a/imulh, hut a minimal improvement (i.e. less than 1 percent) 
was ohserved for all other a/pecIs. 


ri'.o results from the CiMiehition auabsis iiulieateil tliat subtraction of the diffuse component 
leads to very little improvement in the simulation of the radiance by the two ilirecl models evamined, 
with the notable exception of those a/pects perpendicular to the solar a/imuth. In litis case, the 
ci>rrelalion analysis provides only a anirse measure id’ the depreo of association between the direct 
models and the radiance data, l or an explanation of correlation results, it is nece.ssary to examine 
detailed idots of the diffuse and direct data (I’ipure 7). I’or the 2X' and .’b" solar elevation data 
sets, it can be seen that boiii the diffuse and direct radiances increase with slope for those a/pecis 
faciii!’ solar a/imuih. b'or slopes facinp, away from solar a/imuth the direct radiances decrea.se while 
the diffu.se radiances inerea.se. I'he nepalive relationship between the diffu.se and direct radiances and 
the hiph profiortion of the diffuse radiance relative to the direct radiance at bif' a/peet. aecounts 
for the marked increase in the coefficient of determination for the direct radiance models when the 
diffuse component isexiraeled from the global radiance, 

The .second a.s.se.ssment of the significance of the diffu.se component is with reference to band 
ratioing to remove the topographic effect. Band ratioing has been shown to provide a means of 
reducing the topographic effect on remotely sensed data (Vincent, 1^73). As the diffuse component 
of the radiance has been shown to vary with waveletigth. band ratioing of two spectral channels (i.e., 
Channel i/('hannel j) will not le,ul to complete elimination of the topographic induced variations 
(llolben and .luslice IbHO), 


To assess the effect of the diffuse light component on band ratioing of the two channels, the 
ratios were calculated for both the global and direct radiances for all slopes and a/pects (Hgure 8), 
Subtraction of the diffuse component from the radiance data, led to a .''0 percent average decrease 
in the standard deviation in the ratioed values for all a/pecl classes. I’he degree ol reduction was 
greatest, approximately percent, for a/pecIs perpendicular to the principal plane and least, ap- 
proximately .70 percent, for a/peel classes parallel to the principal plane. 
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riu' lai'i’e liciMvaso in iho staiulard deviations of the iliivet lii’ht ratio values rii'in tlie I'toKil 
ladianee ratios, is due to diffuse skylit'lit. wliieli eannot be renuned by simple band ratioinr.. the 
results from this analv sis slu>\v th.it the diffuse eompimeni eontributes siimifieantb to the variation 
in speetral band ratioinp, and I'lvvents ratiiunc. fnun completely removinj'. the topographic effect m 

I the field data. 

I 

•1.5 Assessing the Piffuse I ipht ('ontribution to Simulated 1 aiulsat MSS Radiances. 

^ riie fiekl measurements esamined in this stiub reiuesenl an ideal controlled data set which 

S may help to predict the topographic effect on other sen.sors of multisi'ectral ilata. I aiulsat MSS 

' data are perhai's the most wideh used tv ih' of remotelv sensed multispectral data in earth resources 

I 

' analysis and have been shown to exhibit m.irked topographic effects in the areas ol nipped terrain 

( 

[ (llolben and .lustice. inni) p, effort to assess the importance of the diffuse component on 

' 1 aiulsat radiance ilata, the field measuivil radianees collected in this studv were eonverteil io 1 and- 

I 

‘ sal } (MSS 5 and pixel value equivalents. The simulated pixel values were obtaiiuul by convertinp. 

the haiul-helil radiometer radiances to radiaiu-e values that would be received by the I andsat sensor, 
p.iven specified atmospheric conditions and then iiuantisinp, the radiances accordinp. to the linear re- 
si'on.se of the I aiulsat d analop.ue-to-dipital converter. The method used for this simulation studv 
is detailed in the appendix and was u.sed in a 1 aiulsat simulation study by fucker t fhe I'b- 

jeei of this simuhition studv was to demonstrate the ran.e.e of pixel values li.e., I andsat quanti/ation 
levels), that could be expected from the diffuse component and thereby assess the importance to 
cover classification and niodelinp the topopraphic effect on 1 andsat data, Radiances derived from 
both the RaSO,| and the coarse sand surface were used for this study, flie RaSO,, surface had c. ns 
percent refleetanee which was far hipher than the reflectance from most natural surfaces and in een 
tain cases led to saturation of the simulated satellite sensors, fhe coarse sanil surface had approxi- 
mately oO I'ercent reflectance, fhese two surfaces repre.sented some of the iiiphest radiances from 
natural surfaces that would be obtaiiuui by the satellite for the solar elevations in question and would 
in turn result in the preatest possible topoprai'hic effect. Similarilv’ the raiij’.e of slojU’s examined 
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0 -dO". ivpivsonioil a rargivatcr range (lian is commonly Ibmul even in areas of ruggeil lerrain 
aiul tlierelbre exaggerateil the tiegree of topographic effect that eoulil he i peelcil. To provide 
a more realistic representation of the diffuse liglit topographic effect on l.andsat data, the range of 
simulated pi\el values was calculated for 0 * Mf slopes for each a/pect . as well as for the ct)inplete 
a/pect string ( Table 2), The greatest range in simulated pixel values for the sand surface for the 
0- (H)'’ slopes was for the 0" a/ivet in MSS 7, whereas the maximum range for the ()- slopes 
was } pixel values, 

A maximum range of 3 pixel values for the diffuse light from a (d) percent reflecting surface 
for all a/pects at .slopes of 0 - 30'’ leads to the question of whether consideration of the diffuse 
component should be included in analysis of l.andsat data. The refleetances associated with most 
natural surfaces fall well below oO |iercent and as such the diffuse topographic effeet would be 
even smaller. Given the range of pixel values a.ssociated with any cover type it is unlikely that (he 
diffu.se topographic effeet would be significant for cover type discrimination. However, where 
adjacent surface refleetances ami slopes are high the diffuse component may need to be 
considered. 


.s. DISC USSION Ol- RTSUU'S AND (X)Nn.USIONS 

Tile diffuse light contribution to (he topographic effect was described as a function of solar 
elevation, atmospheric conditions, surface geometry, surface reflectance propeities; and surround- 
ing terrain reflectance. The results presented in this study, although specific for (lie particular 
conditions and surfaces examined, detail a number of general trends and implications common for 
a range of solar elevations (hat can be extrapolated to other studies. Additionally this study 
evaluated the importance of the diffuse light component for three norniali/ation tecliniques. 


Quantification of the diffuse radiances showed a marked topographic effect, tlie greatest 
variation being in the principal plane. Railiances from theBaSOa aiul sand surface were found to 
increase with slope and for tlie rang<’ of surface orientations examined, were strongly influenced 
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by the anisotropic skylitiht ilistrilnition, Railianees rroin liijrli angle slopes were slightly inereased 
by light reriected from siirrouitding terrain. The sand surface used in this study was found to 
exhibit no preferred orientation of seattering. The diffuse radiances cxhibiteil the same pattern 
in both the channels examined, although the diffuse radiances were larger and the topographic 
effect was higher for the photographic infrared channel than for the red ehannel. The higher 
radiances were due to thi wider bandwidth of the infrared channel, Although the diffuse radiances 
varied relatively little compared to the direct radiances, the lu'oportion, of the diffi.se light relative 
to the direct light varieil consiilerably, particidarly for slopes facing away from solar a/imuth. 

Both the diffuse railianees and the diffuse topograidiic effect were found to inerease with solar 
elevation, The diffuse light component under clear sky conditions for the horizontal surface ranged 
from between 10 and 14 percent of the global radiance solar for elevations for 23 to How- 
ever, radiances recorded at low solar elevations were found to have a smaller percentage diffuse 
component for slopes into solar azimuth and greater percentage diffuse component for slojU’s away 
from solar azimuth than at higher solar elevations. 

The contribution of the diffuse light to|iographic effect to the overall topographic effect within 
global radiances, varies with azpect; the diffuse radiances generally increasing with slope for azpects 
facing into sun and decreasing with slope for azpects perpendicular and away from sun. These 
results indicate than an isotropic sky assumption cannot be used to adequately deseribe the diffuse 
component. 

Assessment of the significance of the diffu.se component was undertaken first by examining 
the effect on two direct radiance models. Only slight improvements were found in the correlation 
between the radiance and the theoretical data derived using the Lambertian Model, whereas an 
average improvement of 26 percent in the correlation coefficient was calculated for data derived 
using the non-Lambertian model when the diffuse component was subtracted. 
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|lu‘ signil'icuncc of the ilitTusc component to baiui ratioing was also assessed. Subtraction 
of the iliffuse component Icil to a .‘>0 percent average tiecrease in the standard deviation of 
the ratios associa.ed with the uniform surface, over the range of slopes and a/pecis e\amined. 
Ihc greatest reduction in the standard devaitions was found for slopes perpendicular to solar 
a/iimith. 

I hc fandsat simulation study showed that for clear sky conditions over a range of solar eleva- 
tions associated with typical mid-latitude l.aiulsat passes, that the diffuse light radiances in cases 
of extreme reflectance and slope variation woidd constitute only a maximum 3 pixel variation in 
the I aiulsat measureil global railiances. In this event, it h unlikely that such a maximum variation 
could be taken into consideration to improve cover classification accuracies, by reducing the diffuse 
light toimgraphic effect. It should be understood that areas of shadow will have 100 percent diffuse 
radiance and the higher Landsat uuanti/ativm values associated with these shadowed areas are essen- 
tialb' due to atmos(iheric path radiance. The effect of varying atmospheric path radiance betweeit 
the target and sensor were not examined in this study. 

One methoil used to eliminate the diffuse component and incoming path radiance effects from 
l andsat radiance data is dark object subtraction (Bentley et al.. I ‘>70). I'or this method the lowest 
radiance within the .scene is subtracted from all the radiances prior to ratioing. This method assumes 
an i.sotropic sky distribution and a Lambertian surface and, depending on the range of diffuse radi- 
ance values in the scene, will lead to an over or underestimation of the radiances and subsequent ilis- 
tortion of the resulting ratios, 


Diffuse light cau.ses a topographic effect on remotely sensed data which will vary in significance 
with the application in question. The ground ba.sed study showed that the diffu.se component 
cau.scd a detectable variation in the ratioed spectral data which potentially could be modeled and 
hence the effect eliminated from the data, I'or Landsat studies a maximum possible variation of 
^ pixel values for the conditions specified does not at present warrant further detailed modeling 
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aiul in mosi casos Uic variatHm due to tlio dili'use component will he suhstantiall> less. 
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AIM’I NDIX. CONVI KSION OI 'RADIANa s TO I ANDSA I QUAN ll/AI ION 1 I VI I S 

The liilTuso railiaiK'cs moasurcil by ilu' haiul-licUl lailunuclor wore coiiveiU'il maduMualically to 
radiaiicos lliat woiiKi bo rocoivoil by tlio l.aiulsal imiltispoolral scaiinor (MSS) sonsois. I boso laiii* 
anoos woiv (bon iiuanti/od aooording to iho linoar ivsponso oi llio 1 aiulsat } anaUn'- lo ilipital oon- 
vorUT. 

riio two spootral obamiols of tlto baiul-liokl nulioniolor wore filtoroil to roprosont ttio proposod 
I'homatic Mapper (TM) Iniiuls d and 4 of I aiulsat I) (0,od-0.o‘) jum and 0.7o-(),n0 jimi). Tlioso (wo 
oliannels are sulTioiontly similar to (ho MSS ohannols 5 and 7 (O.o-O.? juin and O.S-1,1 jum) id' I and- 
sal d, to onablo the following simulations. 

Tho Taiulsat d MSS sensors are oalibratod in mw/om' -sr. I'oeonvort (ho hand-hold railiomotor 
radiances from W/m* in (ho I'M bandwidth to tho MSS oalibration, Iho radianoos wore divided In 
10(W/m-’ *■ mw/om’ ), ilividod by tt to oonvorl (i» a storadian moasnro and nuiltipliod b> Iho ratio 
of tho appropriate MSS baiulwiilth to tho I'M bandwidth. This prooodiiro assumes no atmosphorio 
degradation or oontribution to Iho observed signal. Tor this siimilaiion study wo oho.se tho oxtromo 
o.xamplo of (ho maximum MSS rospon.se to tho dilTu.se radianoo emanating from tho surfaoo, from a 
target ol uniform oovor and infinite extent. 

Onoe tho radianoos wore oaloulatod. (hoy wore quan(i/od from 0-127 levels (2' bite) Table 
A provides tho maximum/minimum MSS radianoos nooossarv' lo form tho linoar regression for 
i)uan(i/ing tho radianoo data. 

Table A. Minimum and maximum railianoos timv/om’ - sr) oorrosponding to i.iuan(i/ation levels 0-1 27. 

MSS Ouanti/.ation level 

(’hannol 0 1 27 

4 .04 2.50 

5 .0.) 2.0 

(1 .0.^ l.o.S 

7 .0.) 4.50 

Soma’: lin'l rnsk\ . I'ISO, IVrsoil.il roiiiiluiiiiiMlioil. N,\.SA/li.SI 
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Kv’j’iossioii i-qualu'iis wou’ i-.ili-ulalod for oaoli haiul atul llio sli>|U‘s (m) atui inU'ivcids (b) aiv 
ptvsciUod m lablo H. I'lu' iviirossion ciiualioii lakes ilio lorm of: 
khianti/aljon love' m\ Radiaiu'O i b 

I able It. riie slopes and inleivepts ol' the linear qnanli/alion regression presenled for each MSS 
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SOLAR ELEVATION 29« 



I'ii-.uiv 1. nilTuso RailiaiKV iW'iu’ - sr) (."’o •• MUm) I’lottoil R\ Slope and 
A/'poet - Barium Sulplialo Surrace 




SOLAR ELEVATION - 29° 



SOLAR 1 
NORMAL Y 



Figure 2 , OilTuse Riulianco (W/nr - sr) (.76 - ,‘)0 /um) Plotted By Slope and 
Azpect - Coarse Sand Surfaee 



SOLAR ELEVATION ■ 29° 



SOLAR ^ 
NORMAL ■if' 



180° 


Figure 3a, Reflectance Values for Direct Light (.76 - ,90 juiii) Plotted by Slope and Azpect 

Coarse Sand Surface 



SOLAR ELEVATION » 29“ 



Figure 3b, Reflectance Values for Diffuse Light (.76 - .90 /am) Plotted by Slope 
and Azpect - Coarse Sand Surface 



SOLAR ELEVATION - 23° 



I'igure 5a. Diffuse Radiance (W/nr - sr’ '.63 - ,6‘) fim) Plotted by Slope and 
Azpect - Barium Sulphate Surface 






SOLAR ELEVATION - 39® 



Mtunv 5c. DilTuso Riuliaiice (W/nr - sr) (.63 - .66 /itm) Plotted by Slope and 
A/,pect - Barium Sulphate Surface 




Figure 5il. DilTuse Radiance (W/m“ ■ sr) (,76 - ,90 /am) Plotted by Slope and 
A/pect - Barium Svilpliate Surface 



SOLAR ELEVATION • 39® 



I’igiire (nt, Percentage DilTuse Sunlight (.63 - .6‘) Mm) Plotted by Slope and A/.pect 

Barium Sulphate Surface 





fR/RED RADIANCE 
RATIO VALUE 



iMguiv S, RailiaiKv Values aiul Slaiulard Doviations ('alculalod 

iVoiu Cilobal and nirocl Hadiaiu’os vs. .V/iu'ct 



